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ABSTRACT 

A pilot project to establish an operational phase stable very long baseline interferometer (VLBI) for 
geophysical studies is described. The technique of phase transfer between distant oscillators has 
been demonstrated in an experiment using the Hermes satellite. Based on the results obtained, a 
phase link employing two-way transmissions of multiple pilot tones via standard satellite commu- 
nications channels is designed. The project will employ the Anik-B satellite to provide a phase link 
between ARO in Ontario, Dominion Radio Astrophysical Observatory (DRAO) in Penticton, and 
Naval Research Laboratory (NRL) in Maryland Point. VLBI observations of a selection of strong 
compact radio sources will be carried out at two frequencies using the Canadian analog VLBI re- 
cording system. By correcting the observed fringe phase using the satellite link phase, a major source 
of uncertainty in VLBI geophysical measurements can be removed. 

The instrument will be used to monitor Earth’s rotation, the nutation of the Earth’s axis, and other 
geodynamical observations. One objective of the project is to evaluate the feasibility of a dedicated 
phase coherent geophysical VLBI for UT1 monitoring to replace the PZT’s. 
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INTRODUCTION 

The possibility of phase coherent long baseline interferometry was pointed out by Yenetal. [1977] 
following a series of experiments in which the Hermes communications satellite was used to provide 
a data link between antennas in Canada and the United States to operate a long baseline radio inter- 
ferometer in “real time”. Following the success of this effort, an application by J. L. Yen and W.H. 
Cannon to the Federal Department of Communications, Ottawa, Canada, was approved for the use 
of the Anik-B communications satellite for purposes of constructing a phase coherent long baseline 
interferometer for geophysical applications. The Anik-B communications satellite was launched on 
December 9, 1978, and is owned and operated by Telesat Canada. The satellite’s 12/14 GHz trans- 
ponder capacity has been leased from Telesat Canada by the Federal Department of Communica- 
tions and has been made available, free of charge, by the Department of Communications to 
qualified users. 

The initial objectives of this research program are to investigate the potential of domestic telecom- 
munications satellites for conducting phase coherent long baseline interferometry and to test the 
capability of a phase coherent long baseline interferometer for monitoring short term variations in 
the Earth’s rotation rate and other geophysical studies. 

In carrying out these studies, the Anik-B satellite will be used to establish a phase link between the 
following observatories: 

1 . ARO, Algonquin Park, Ontario, Canada 

2. DRAO, Penticton, B. C., Canada 

3. NRL, Maryland Point, Maryland, U.S.A. 

The astronomical receivers at each observatory will be dual frequency, 1 1 cm and 21 cm, systems 
using cavity backed spiral feeds. The astronomical data will be recorded using the Canadian analog 
format on video cassette tape recorders with a 3.5 MHz instantaneous bandwidth at -3dB. The long 
baseline interferometry recording terminals used at each observatory as well as the correlator facility 
used to obtain the interferometer fringes will be provided by the Radio Astronomy Section of the 
Herzberg Institute for Astrophysics, Ottawa, Canada. Some general characteristics of this interfer- 
ometer are summarized below in tables 1 and 2. 


Table 1 


SITE 

ASTRONOMICAL 

ANTENNA 

DIAMETER 

PHASE LINK 
ANTENNA 
DIAMETER 

PHASE LINK 
TRANSMITTER 

ARO 

46 meters 

10 meters 

20 watt TWTA 

DRAO 

26 meters 

4.5 meters 

20 watt TWTA 

NRL 

26 meters 

26 meters 

20 watt TWTA 
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Table 2 


INTERFEROMETER 

PARAMETER 




DRAO-ARO 

DRAO-NRL 

ARO-NRL 

X 

2970 tan 

3160 km 

190 tan 

Y 

-720 tan 

-1260 km 

-540 km 

Z 

-250 tan 

-880 tan 

-620 km 

TOTAL LENGTH 

3060 km 

3510 km 

840 km 

EQUATORIAL LENGTH 

3050 km 

3400 km 

570 km 

EAST LONGITUDE 

-13° 

0 

OJ 

(N 

1 

1 

H 

0 

DECLINATION 

-5°. 

0 

f— 1 
1 

o 

00 

1 

MIN. FRINGE 
SPACING 11 cm 

7V4 x 10~ 3 

6 '.'4 x lo" 3 

2 7 V 0 x 10 -3 

MIN. FRINGE 
SPACING 21 cm 

14 '.'1 x 10 -3 

12'. ‘2 x 10“ 3 

5 1 5 x 10 -3 


The entries in table 2 have been rounded off and are therefore only approximate values. 


THE OPERATION OF A SATELLITE PHASE LINK 
General Description 

Radio interferometry ideally requires local oscillator signals at each antenna site whose phases re- 
main coherent over extended, ideally infinite, intervals of time. In connected element interferom- 
etry, this property is easily achieved by distributing local oscillator signals via cable or microwave 
links from a central source. In long baseline interferometry (LBI) the separation between the an- 
tennas generally precludes the generation of phase coherent local oscillator signals by such means, 
and local oscillator phase coherence at each antenna site is traditionally achieved by phase locking 
the local oscillator, independently at each site, to separate frequency standards whose absolute 
phase stability is known to be very high. 

Since the phase observable is extremely sensitive to geodetic and geophysical phenomena consider- 
able effort has been devoted to improving the intrinsic (i.e., uncorrupted by atmospheric propaga- 
tion effects) phase stability of long baseline interferometers. Part of this effort has been directed 
towards the development of frequency standards, particularly hydrogen masers, of extraordinarily 
high absolute phase stability [Vessot et al., 1977] . However, such devices are becoming both ex- 
pensive and bulky and will require operation in very carefully controlled environments to achieve 
their maximum potential [Vessot et al., 1977] . 
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While absolute phase stability over extended periods is essential for many applications such as the 
Doppler tracking of spacecraft, it is inessential for long baseline interferometry applications which 
carry with them the somewhat weaker requirement of phase coherence or relative phase stability 
over extended periods. The satellite phase link is a technical development which is intended to ex- 
ploit the opportunity presented by the somewhat relaxed requirements of the interferometer appli- 
cations. There appear to be no fundamental physical impediments preventing satellite phase links 
using inexpensive and robust rubidium frequency standards or even high quality crystal oscillators 
from achieving long baseline interferometer operation with intrinsic phase stability comparable to 
that achievable using super stable hydrogen masers. 

The operation of the Anik-B satellite phase link is illustrated in figures 1 and 2. Figure 1 illustrates 
the recording procedure; figure 2, the playback procedure. In the Anik-B experimental program, 
the recording and playback operations will remain distinct. However, in an operational system, they 
would be combined by sending the recorded signal from one antenna to the other, whereupon the 
mixing shown in figure 2 would be performed, thereby achieving real-time operation of the phase 
link. 

In the present configuration of the Anik-B phase link, the phase of the astronomical local oscillators 
(p 1 and $2 at antenna sites #1 and #2, respectively, are used, along with four phase-locked local 



Figure 1 . 
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4-2 - *1 


Figure 2. 


oscillators to generate two pairs of signals. The high (upper) frequency pair, denoted by a super- 
script U, and low (lower) frequency pair, denoted by a superscript L, are separated by 60 MHz in 
each case. These signals with phase K^ 1 <t> x , (p l at site #1 and K 2 <t > 2 , K 2 <A 2 at site #2, where 
, Kjf and K 2 , K 2 are known rational numbers with > K\ and K 2 > KL , are transmitted to 
Anik-B via the 14 GHz up-link channel. The phases of the signals arriving at Anik-B from site #1 

are K} 1 [<p l - </> G1 - 4>ui ] > [01 “ <t> gi - 0mi 1 and from site # 2 are K 2 [<p 2 - <p C2 ~ <P U 2 ] > 

K 2 [</> 2 — 0 G2 — 2 ] , where 0 G1 , </> G2 refer to the geometric phase delays and 0 M1 , 0 M2 refer to 

the propagation medium phase delays at the local oscillator frequency between Anik-B and the 
antenna sites #1 and #2, respectively. 


In transitting the satellite, the translation oscillator on board Anik-B converts the up-link signal 
at 14 GHz to a down-link signal at 12 GHz, subtracting the phase of the onboard translation oscil- 
lator <p T . The phases $ 21 , 0 21 of the signal arriving at site #1 from site #2 are therefore: 

021 = K 2 [02 — 0G2 — 0M2 1 ~~ ^2 I0G1 + 0M1 ^ ~ 0T 
021 = ^2 [02 — 0G2 — 1 M 2 ] — K 2 [0 G2 + 0M1 1 — 0T > 
and the phases 0 ^ 2 , <p^ 2 of the signal arriving at site #2 from site #1 are therefore: 

012 = ^ [ 01 — 0G1 _ 0M1 J ~~ ^2 [0G2 + 0M2 1 _ 0T 
012 = ^1 [01 _ 0G1 _ 0M 1 1 _ ^ [0G2 + 0M2 1 ~~ 0T > 
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where the factors , Kj and are introduced to account for the fact that the up-link and 

down-link frequencies are different. 

The three stages of mixing which occur at each ground station as shown in figure 1 produce record- 
ed signals with phases (hj and <J> 2 at stations #1 and #2, respectively, where 

= 1*2 -*2 11*2 -<P G2 -[K 2 u '-K 2 L '][0 G1 + 0 M1 ] - [Kj 

$2 = [Ky - K\ ] [0J - 0 G1 - 0 M1 ] - tKy' - K\') [0 G2 + 0 M2 1 - [K« -K[]0 2 . 

Now Ky -K\ = Ky - and Ky — = Ky — , since the frequency separation on the up- 

link and down-link are identical. Hence, the recorded phases at each antenna site are: 

4>J = [Ky - ] [0 2 - 0J - 0 G - 0 M ] 

$2 = _ K 1 1 f^l — *^2 ^ 


where 


0 G 

^G1 

+ 0G2 


= 0M1 

+ ^M2 • 


Following this, it is clear that the playback procedure shown in figure 2 produces the local oscillator 
phase 0 2 — 0j . 


Timing Requirements 

The timing requirements for the data manipulation in a satellite phase link of the Anik-B type can 
be estimated from the expressions for the recorded signals and<f> 2 . The recorded signals and 
<J >2 are seen to depend on the sum of the phase difference 0 2 — 0 t and the phases 0 G and 0 M all of 
which, in this treatment, have been referenced to the local oscillator frequency co Q . 

For local oscillators phase locked to even relatively poor frequency standards with frequency insta- 
bilities as large as 10 -1 ^, the phase rate |-^ (0 2 — 0j )| is of the order of 10 -1 ® co D radians sec -1 . 
For an observing frequency of 3 GHz (co 0 = 6n x 10 9 ), it follows that (0 2 - 0j)| = 2 radians 
sec -1 and a timing error of e seconds introduces an error in the measurement of the order of 2e 
radians. 

The phase rate is a Doppler shift due to the motion of the satellite relative to the ground sta- 

tions, changing the total geometric path length of the phase link. Yen [1977] has shown that for a 
satellite in a geosynchronous orbit inclined at 0.2° to the equator that 
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and that 


0.2586 sec < ^2. 

co n 


< 0.2588 sec 


1 


co 


o 


dt 


| < 5 x 10" 9 


It follows that the phase rate |d0 G /dt| is the order of 5 x 10" 9 co Q radians sec -1 . For an observing 
frequency of 3 GHz (co 0 = 6n x 1 0 9 ), Jd0 G /dtj = 10 2 radians sec -1 and a timing of e seconds in- 
troduces an error in the measurement of the order of 10 2 e radians. 


The phase rate |d0 M /dt | is due to temporal variations in the propagation medium along the path of 
the phase link and will be due to temporal variations in the tropospheric refractive index and the 
ionospheric integrated electron density. Reasonable upper limits can be placed on d$ M /dt by con- 
sidering the magnitude of the tropospheric and ionospheric phase delays and the typical time scales 
of their variations. 


<p M can be decomposed into a tropospheric contribution $ MT and an ionospheric contribution 0 MI , 

= ^MT + ^MI ■ 

If L t is the excess path length over vacuuo imposed on the phase link by the presence of the tropo- 
sphere (the troposphere imposes roughly 2 meters of phase delay in the zenith direction over each 
antenna site) then 

I^mt I = Sfo |dLr | radians sec -i 
dt c dt 

Since the phase link antennas pointing to Anik-B have reasonably large zenith angles, we will take 
Lj ~ 8 meters if we assume that Lj can at most change by as much as 10 percent (80 cm) in an 
interval of the order of 1000 seconds (15 minutes), then 

I 0MT| ^ 2.6 x 10“ 12 to_ radians sec -1 . 
dt 

For an observing frequency of 3 GHz (co Q = 6ir x 10 9 ), |d0 MT /dt| < 5 x 10 -2 radians sec -1 and a 
timing error of e seconds introduces an error in the measurement of the order of 5 x 10~ 2 e radians. 


The ionospheric contribution <p m is given [Rogers and Evans, 1968] by 


^MI 


2 

C0 o 


2i re 2 


m e c 


f L N e d«, 


where e is the charge on the electron (esu), m e is the mass of the electron, and c is the velocity of 
light; N e is the number of density of free electrons and the integral is taken along the path of the 
phase link. Ionospheric studies indicate that 0 M j is generally in the range 
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10 11 ^ _ 3-5 x 10 12 

> <t> u 1 > radians, 


CO 




in the zenith direction at mid-latitudes. Since the phase link antennas pointing to Anik-B have 
reasonably large zenith angles, we will take as a worst case 


4x 10 11 


OOr 


> *MI > “ 


15 x 10 12 


radians. 


COc 


The ionosphere is spatially and temporally quite variable with travelling ionospheric disturbances 
(TID’s) in the F region and sporadic E disturbances in the E region modulating the integrated elec- 
tron content by as much as ±2 percent in an interval as short as 10 minutes. In addition, most 
dramatically at sunrise, but also at sunset, the integrated electron content changes by as much as an 
order of magnitude in an interval as short as 3 hours. 


It follows from this that 


,^MI I 

dt 


1-5 x 10 9 


< 


co„ 


radians sec ' . 


For an observing frequency of 3 GHz (co 0 = 6n x 10 9 ), |d0 MI /dt| <8x 10 -2 radians sec -1 and a 
timing error of e seconds introduces an error in the measurement of the order of 8 x 10 -2 e radians. 


From the above analysis, it is apparent that synchronizing the data on playback to an accuracy of 
e = ±10 -4 seconds is sufficient to insure that errors due to timing do not contribute more than 10 -2 
radians ( — 0?5) of phase error to the measurement. 


The phase rate d<p G (t)/dt is caused by the motion of the satellite at time t, where t refers to time 
kept by a set of synchronized clocks at rest in the rotating frame of the Earth. Phase link signals 
transitting Anik-B at time t arrive at station #1 at time t + d(t) and at station #2 at time t + d(t) + 
5(t), where 5(t) represents the inequality of the travel times from Anik-B to ground stations #1 and 
#2, respectively. Because of retarded electro-magnetic effects, the Doppler shift phase recorded at 
station #1 <p G is given by 


4(t) 



dLc(t' -d(t')) dt , 


and the Doppler shift phase recorded at station #2 <p G is given by 




dLg(t' -d(t') -5(t')) 
dt 


dt' , 


where t 0 is an arbitrary epoch and L G (t) is the geometric length of the phase link path. 
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Clearly, 

0 q ( t ) =£ 0^(0 ; 

however, 

^(t) = 4 >q (t + 5(t)) . 

So, to cancel the effects of the satellite Doppler shift in the playback procedure, it is necessary to 
delay the signal from tape recorder #1 relative to the signal from tape recorder #2 by an amount 
8(t). For simplicity, this delaying operation was omitted from figure 2. 

For ground stations in the beam of Anik-B, 5(t) is about 6 msec. However, we see from the analy- 
sis above that 5(t) need only be known with an accuracy of the order of ±0.1 msec. This implies 
that the location of Anik-B relative to the ground stations must be known with an accuracy of only 
±30 km. 


Phase Link Performance Possibilities 

A discussion of the potential phase stability of the Anik-B phase link can be partitioned into two 
regimes: the first pertaining to time scales of roughly 10 2 seconds or less over which environmental, 
atmospheric, and ionospheric conditions are relatively constant; and the second pertaining to time 
scales of roughly 10 3 seconds or more over which environmental, atmospheric, and ionospheric 
conditions are changeable. 

On short-time scales, the phase stability of the Anik-B phase link will be largely dependent on sig- 
nal to noise ratios within the link itself and on the general stability of the phase lock loops in the 
electronics. 

On long-time scales, the phase stability of the Anik-B phase link will be limited by environmental 
factors which, if left uncontrolled, will considerably degrade the long-term performance of the 
phase link. In the case of Anik-B, there are two principal limitations which will contribute to the 
degradation of the link’s performance. 

The first of these will be the largely temperature-dependent instrumental phase shifts occurring 
principally in the signal cables of the ground stations. These phase shifts will in general not be equal 
at both ground stations and so will not be cancelled by the overall operation of the phase link. 
Fortunately, a variety of techniques are available to either stabilize, eliminate, or calibrate these 
phase shifts, and it would seem that a persistent and systematic attack on this problem could even- 
tually reduce these effects to something of the order of ±5°-10° of phase rms. 

The second long-term limitation on the performance of. the Anik-B phase link is imposed by the 
ionosphere. The ionosphere degrades the performance of the phase link as the result of a small 
second order dispersive effect arising from the dual frequency nature of the two-way link. A dis- 
persive phase shift developing between the two transmitted frequencies traversing the ionosphere on 
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the up-link at 14 GHz is in general not exactly cancelled by a corresponding dispersive phase shift 
between the two received frequencies transversing the ionosphere on the down-link at 12 GHz. The 
small non-cancelling dispersive phase shift denotes 50 MI has been shown [Cannon, 1979] to be 
given by 


where 




a K u - K l 
1-q ’ K u K l 


A 1 (t)-A 2 (t) 




] , 


K u = 


k l = 


+K 


L 

2 


2 


K u + K l u> t 1 /ccxV 

k u k l ZT q k u k l 


A^t) 

A 2 (t) 


27re 2 
m e c Ll 


N e d£ 


27re 2 
/ 

m e c L 2 


N e d£ . 


In the above formulas cu T is the frequency of the onboard translation oscillator or equivalently the 
frequency difference between the up-link and down-link signals. The paths Lj , L 2 for the line 
integrals defining quantities Aj(t), A 2 (2) are taken along the phase link from the ground station to 
the satellite over stations #1 and #2, respectively. 

In general, the quantity [Aj (t) — A 2 (t) ] will be nonvanishing and will usually attain its largest 
values when the interferometer baseline straddles the day/night terminator on Earth. In this worst 
case situation, an order of magnitude difference may easily appear between Aj(t) and A 2 (t). For 
typical values of Aj(t), A,(t), the second order ionospheric effect could produce a dispersive phase 
error 8(p m as large as ±ir for several hours as the terminator passed through. However, aside from 
the day/night transition periods the quantities Aj(t), A 2 (t) are more likely to differ by factors of 2 
to 3 or less. Under these circumstances, S0 MI can be expected to be typically of the order of ±7t/4 
rms. In addition, under both these circumstances, some portion of this second-order ionospheric 
effect is predictable on the basis of standard ionospheric models and may be removed in the data 
processing stage. 

Finally, although the Anik-B program has no provision for it, it is of course possible to calibrate 
these dispersive ionospheric phase errors by the addition of a third frequency to the two-way phase 
link. With the bandwidth available (200 MHz) on Anik-B the dispersive ionospheric phase errors 
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could be calibrated with an accuracy of only ±10 percent. However, this would reduce typical 
values of 8<t> m to something of the order of ±5° rms. 

These results are summarized in figure 3 which shows, for comparison, the operational and theoreti- 
cal performance of the VLG-1 1 hydrogen maser [Vessot et al., 1977] and the operational perform- 
ance of the Hermes (CTS) satellite phase link [Waltman, Knowles, et al., 1979] . It would appear 
that the performance of the Hermes (CTS) satellite phase link which degrades rather badly between 
time scales of 10 2 seconds and 10 3 seconds and longer is perhaps due to uncalibrated instrumental 
and environmental effects. If these are adequately remedied, a satellite phase link performance 
should be able to follow trends given by the solid lines in figure 3 , achieving an equivalent frequency 
stability of the order of parts in 10 16 on time scales of the order of 12 hours, and parts in 10 17 on 
time scales of the order of 4 to 5 days. 


APPLICATION OF THE ANIK-B PHASE LINK TO THE MEASUREMENT OF VARIATIONS 
IN EARTH ROTATION 

General Method of Measurement 

It is intended initially to apply the Anik-B phase link interferometer to the special purpose of mea- 
suring short-term variation in the Earth’s rotation rate. This is a good application of a satellite 


12 hr 1 wk 

1 sec 1 min 1 hr 24 hr 1 mon 1 yr 10 yr 



10° 10 1 10 2 10 3 10 4 10 5 10 6 10 7 10 8 10 9 
AVERAGING INTERVAL, sec 


Figure 3. 
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phase link interferometer as such measurements depend strongly on the longer-term phase stability 
of the interferometer rather than on its short-term phase stability [Broten, 1979] . 

The Earth’s variable sidereal rotation rate £2(t) is related to its mean sidereal rotation rate 52 by the 
small dimensionless parameter m 3 (t) 

fi(t) = fi 0 [l+m 3 (t)]. 

The quantity most readily measured by observing successive times of meridian transits of stars 
against an atomic clock is r(t), the amount by which the variably rotating Earth lags behind a hypo- 
thetical uniformly rotating Earth. 

r(t) = - / m 3 (t') dt' . 

At any time t, the phase of the interferometer can be represented by <f>(t) where 

4>(t) = 0 g (t) + 0 c (t ) + <p A1 (t) + 0 A2 (t) + 0j(t), 

and where 

• 0 c (t) is the so-called clock phase arising as a result of relative phase shifts between two 
local oscillator signals 

• 0 A1 (t), 0 A2 .(t) are t ^ ie resu lt °f relative phase delays on the astronomical signal due to the 
presence of the troposphere and ionosphere respectively 

• 0j(t) is the so-called instrumental phase arising as a result of uncompensated phase delays 
between the receiving/recording systems at each interferometer terminus 

• 0 g (t) is the geometric phase. 

For an interferometer of baseline length b, east longitude A b , declination 0^ observing at frequency 
oj q , a source with right ascension a s and declination 5 $ at Greenwich sidereal time 6, the value of 
the geometric phase is given by 

co b , 

0 g = — ~ [sin0 b ' sin 5 S + cos 0 b cos 5 S cos(0 + A b - a $ )] . 

Beginning at an epoch t 0 and carrying a continuous “fringe count” forward in time for an interval 
t 0 < t < t 0 + r yields the observable function of time 


[4>(t) - <J>(t 0 )] = 0 g (t) - 0 g (t o ) + 0 c (t) - 0 c (t Q ) + 0^(1) - 0 A1 (t o ) + 0 A2 (t) - 0 A2 (t o ) + 0,( t) - 0,(t o ). 
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Repeating the measurement, an interval of n sidereal days later for a corresponding interval t Q + nT 
< t < t Q + nT + r yields 

[4>(t + nT) — <h(t 0 + nT)] = 0 g (t + nT) - 0 g (t o + nT) + 0 c (t + nT) - 0 c (t Q + nT) 

+ <?>A1 (t + nT) - 4> m (t Q + nT) + ^(t + nT) - + nT) 

+ 0j(t + nT) - ^j(t 0 + nT) , 

where the interval T represents the a priori estimate of the length of the sidereal day. Differencing 
these two data sets gives 

{ [<f>(t + nT) — <J>(t Q + nT)] — [4>(t) — 4>(t 0 )] } - { [0 c (t + nT) - 0 c (t Q + nT)] 

- [0 c (t) -0 c (t o )]j - jt^Ct + nT) -0 A1 (t o +nT)] - [0 A1 (t) -0 A1 (t o )]| 

- | [^(t + nT) - ^(to + nT)] - [0^(0 - ^Oq)] j - | [0 x (t + nT) 

- <t>i( t G + nT)] - [ 0j(t) - ^(t 0 )] | = 0 g (t + nT) - <j> g (t) - K n 
where the quantity K n is a constant for fixed n given by 

K n = 0 g (t o + nT) - 0 g (t o ) . 

Each term within the braces consists of the difference of two terms within brackets and the brack- 
eted terms are all observable quantities within some experimental error. 

• [4>(t) — 4*(t 0 )] can be measured by counting fringes at the output of the correlator 

• [<£ c ( t) — 0 c (t Q )] can be measured by the Anik-B phase link 

• [0 A1 (t) — 0Al(t Q )] can be estimated using ground based pressure temperature and humid- 
ity measurements along with a standard atmospheric model 

• [0A2^ 0a2^o^ can nrieasured by dual-frequency, 11 cm and 21 cm, astronomical 

observations 

• [0j(t) — 0j(t o )] can be measured by the use of an instrumental phase calibrator system. 

It can be shown that to first order in small quantities: 
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0 g (t + nT) - 0 g (t) = | — A0(t, nT) cos 0 b ' cos 5 g sin (0 + A b - a g ) + A0 b '(t, nT) 

[cos 0^ sin 5 g — sin0 b ' cos 5 S cos (0 + A b — a s )] + A5 g (t, nT) 

[sin 0/ cos 5 g — cos 0^ sin 6 g cos (0 + A b — a g )] - AA b (t, nT) 
cos <j>J cos 5 g sin (0 + X b — « s ) + Aa g (t, nT) cos 0 b ' cos 5 g 
sin (0 + X b — a & ) + Ab (t, nT) -g [sin 0^ sin S g + cos 0^ cos 8 $ 
cos (0 + X fe -a s )]| , 

where Ab(t, nT), A0 b '(t, nT), AX b (t, nT) represent increments in the interval t < t' < t + nT to the 
baseline scalar length, declination, and east longitude respectively ; and where Aa s (t, nT), A5 S (t, nT) 
represent increments in the interval t < t' < t + nT to the source coordinates of right ascension and 
declination respectively. The quantity A0(t, nT) is defined in such a way that the total increment 
to the sidereal angle of the baseline which occurs in the interval t < t' < t + nT is 2nn + A0(t, nT). 
Hence, 


A0(t, nT) = - Ar (t, nT) ; 


or 


A0(t,nT) = n o / +nT m 3 (t')dt' . 

t 

Table 3 summarizes the properties of the partial derivatives of the differenced interferometer phase 
with respect to the variations in the above parameters. 

For measurement of variations A0 in Earth rotation, it will be useful to adopt an observing strategy 
which conce ntrates on low (5 $ < 10°) dec lination sour ces. Since 0 b ' = -5° on the baseline 
DRAO-ARO and 0 b ' = -14° on the baseline DRAO-NRL, it is clear that by observing low declina- 
tion sources the interferometer can be rendered relatively insensitive to the variations A0 b ' and A5 g . 
The variations A0 b ' arise as a result of polar motion parallel to the meridian of the interferometer 
and the variations AS g arise as a result of precession and nutation in the interval t < t' < t + nT. 
Furthermore, since polar motion mj(t) m 3 (t) increments the interferometer baseline longitude by 
AA b where 

sin 0v,' 

AA b = — [mj(t) sinX b - m 2 (t) cosA b ] , 

COS 

we see that these low declination baselines are also rather insensitive to variations AA b caused by 
polar motion. The values of the parameters A0 b ', AA b , A5 $ to which the interferometer is relatively 
insensitive but which are required for data processing can be obtained from sources such as the 
Bureau International de THeure (BIH) and the standard IAU theory of the precession and nutation. 
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Table 3 


ASTRONOMICAL 

PARTIAL DERIVATIVE 


GEOPHYSICAL 

CONSTANT TERM 

w b 

(0 b 

PARAMETER 

VARIATION 

c 

o 

X 

c 

AMPLITUDE 

SIDEREAL TIME 
DEPENDENCE 

A0 

— 

- cos cos & s 

sin 0 

A*f> ^ 

cos d> ' sin 5_ 
b s 

- sin cos 5 

■O s 

cos 0 

A6 

s 

sin d>' cos <5 
b s 

— cos it' sin 6 
b s 

cos 0 

AX b 

— 

- cos <b' cos 6 
b s 

sin 0 

Act 

s 

— 

cos <b' cos 6 
b s 

sin 0 

Ab 

— sin 6' sin 6 
b b s 

ir cos cos 6 

b b s 

cos 0 


Of the three remaining parameters to which this data is sensitive, we see that A 8 and Aa s are insep- 
arable. It will therefore be necessary to use existing IAU models for the precession and nutation to 
theoretically compute the changes in right ascension Aa s of the source position and finally to solve 
for quantities A0(t, nT) and Ab(t, nT) separating them on the basis of their different sidereal time 
dependences. 

The entire measurement procedure described above is insensitive to errors in the adopted a priori 
values for the parameters of the interferometer. The effect of such errors on the measurement is 
made to vanish by differencing the data from day to day. 


Performance Possibilities for Measurement of Earth Rotation 

The rms measurement errors y/ 86 2 andVSb 2 in the parameters A 6 and Ab are dominated on the 
time scale of days to months by the rms error V 8<p 2 in the measurement of the differenced inter- 
ferometer phase rather than by systematic errors in the modelling of the quantities A <^' AA^ Aa $ 
A5 S . The latter quantities, particularly Aa g , A8 S , will become important sources of error on time 
scales of decades to centuries as errors in the models for the precession and nutation will inevitably 
propagate into errors in the measurement of earth rotation. 

The parameters Ad and Ab are quite well separable on the basis of 12 hours or more of differenced 
data; so, for data sets of 12 hours or more, we can neglect the covariances between Ad and Ab and 
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obtain a first order estimate of their respective measurement errors directly from their partial deriv- 
atives. We can take 


or equivalently 


and 


VSF a y/2 ~TT 

"o b 


y/br 2 = y/2 


n o w o b 



y/tt? = y/2 — ^ , 

where the factor y/2 is to approximately allow f or th e fact that the value of the partial derivative 
varies sinusoidally with sidereal time and where y/ 8 / 1 , given by 



\/2(50 2 + 50 2 + 8<t> 2 Al + 8<j> 2 A2 + 50 2 ) 


is the instantaneous (one sample) rms measurement error of the differenced interferometer phase. 
The factor y/2 in the above formula is required to account for the fact that we are dealing with dif- 
ferenced data. 


We shall assume that the stochastically varying phase errors in the measurements of [<f>(t) — <f>(t 0 )] , 
[0 c (t) - 0 c (t Q > ] , [0^ (0 “ ^Al(t 0 )], ~ ^A2^o^> and I^I^) “ 0i(t o )]> having variances 

502, 502, 502 j , 502 and 502 , respectively will be characterized by time scales of independent 
fluctuations T o’ T c> T A1 ’ T A2’ and Tj, respectively. It foll ows t hat a least squares adjustment to a 
data set of length r will be corrupted by an rms phase error y/s^ir) given approximately by 

\/W(T) = y /- T (T q H 2 + T c 50j +T A1 5^ + TA2 + Tj 50* '. 


v /50 2 ' is the rms measurement error on the phase of the correlation function. For strong sources 
vV' is generally quite small, of the order of ± 10°-15° and since each integration provides an inde- 
pendent sample of this random variable the time scale T 0 is of the order of seconds. 

M is the rms measurement error on the clock phase and wil l be dominated on the time scales 
under consideration here by instrumental phase delay variations 50 2 j in the cabli ng an d electronics 
of the ground stations and by the dispersive ionospheric phase delay variations 50 2 2 which affect 
the phase link operation. 


T Me = 


T cl + T c2 6 * 2 c2 ■ 



ADHERENCE AND COHERENCE IN NETWORKS AND PLANS 


477 


We have seen that could be as large as 7r when the day/night terminator passes through but 

could be expected under other conditions to be typically of the order of ir/4 with a time scale of 
fluctuations T c 2 of the order of 30 minutes, the coherence time of the short-term temporal varia- 
tions of the ionosphere. 


In the absence of phase calibrators, the repetitive component of the diurnal variations in the 
instrumental phase delays at each ground station will be cancelled by the phase differencing 
procedure. The residual variations in the instrumental phase delays can be expected to be of 
the order of 2 nsec rms with a time scale of fluctuation of several hours. For a frequency of 
3 GHz, we can take the uncalibrated rms variations to be of the order of \2ir with a 

time scale of fluctuations of the order of 2 hours. 


The addition of phase calibrator systems to the Anik-B phase link could be expected to greatly 
reduce this source of error, particularly since the operation of the phase link merely requires the sta- 
bilization of the instrumental phase delays; a much weak er req uirement than their actual measure- 
ment. It would seem that phase calibrators could reduce >/ 50§ 1 1:0 quite small values with something 
of the order of tt/10 reasonably attainable. In such a case, however, T c j would remain the same, of 
the order of 2 hours. 


To first order, we can apply these remarks to the error source which arises from phase delay 

variations in the cabling and electronics of the astronomical receiver/recorder syst em. T hus, we can, 
for purposes of our estimate, take = 1 2n without phase calibrators and y / 502 s ir/lOwith 
phase calibrators and Tj = 2 hours in each case. 




„ A1 is the rms phase error arising from the failure to properly correct for the phase delay effects 
of the troposphere on the astronomical signal. Using ground based pressure, temperature, and 
humidity data along with a standard atmospheric model it is possible to compute the tropospheric 
phase delay with an accuracy of 5 cm rms at the zenith and 10 cm rms at a zenith angle of 60°. The 
additional use of dual frequency water vapor radiometers can apparently reduce these errors to 
something of the order of 2 cm rms at the zenith and 4 cm rms at a zenith angle of 60° . Since a 
long baseline interferometer almost always has one or other of its antennas pointing at a large zenith 
angle, we shall take the 60° zenith angle figures as being typical. For a frequenc y of 3 GHz, we can 
assume \J§(p 2 M ^ 27r without the water vapor radiometer measurements and \/d(f > 2 A1 ' — 0.87T with 

the water vapor radiometer measurements. The time scale of fluctuations of these quantities will be 
taken to be 30 minutes, roughly the coherence time of the short-term variations of the atmosphere. 


y/t 


5 0 A2 is the rms phase error arising from the failure to properly correct for the phase delay effects 
of the ionosphere on the astronomical signal. Dual frequency observations at 1 1 cm and 21 cm 
allow removal of the effects of the variations in the differential ionospheric phase delay along the 
beams of the antennas with an rms error of roughly 1 percent. For a 3 GHz signal, the magnitude 
of the ionospheric phase delay is of the order of 2 meters at zenith angles of 60°. Although large 
scale excursions in differential integrated electron content occur near sunrise and sunset, the differ- 
ential integrated electron content can be expected to be typically much smaller. We shall assume a 
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25 percent differential integrated electron content between the two antennas as being typ ical. If 
we again assume a 60° zenith angle as being typical of the interferometer operation, then-y/g^^ 

0.1 7r at 3 GHz. The time scale of independent fluctuations of this quantity can be taken to be of 
the order of 30 minutes, the coherence time of the short term ionospheric variations. 

The results of this error budget are summarized in table 4 along with the quantity >/5(^ 2 (r) which 
results, r is measured in hours of data fitted. 

Figures 4 and 5 based on the error budget of table 4 indicate that the addition of phase calibrator 
system to the astronomical receiver/recorders and to the phase link ground stations will greatly im- 
prove the performance of the interferometer. However, beyond this, the further addition of water 
vapor radiometers and a third frequency in the phase link to remove ionospheric effects can offer 
only marginal improvement by perhaps a further factor of two. 


Table 4 


ERROR 

SOURCE 

TIME SCALE 
OF 

INDEPENDENT 

FLUCTUATIONS 

APPROXIMATE RMS AMPLITUDE 

DUAL 

FREQUENCY 

OPERATION 

PHASE 

CALIBRATORS 

ADDED 

WATER 

VAPOR 

RADIOMETERS 

ADDED 

PHASE LINK 
THIRD 
FREQUENCY 
ADDED 

Sf> 2 

10 sec. 

0.08 TT 

0.08 TT 

0. 08 tt 

0. 08 TT 

5 *cl 

2 hr. 

12 7 r 

0. 1 TT 

0.1 TT 

0.1 7 T 

5 *c2 

2 hr. 

0.25 TT 

0.25 tt 

0.25 tt 

0.02 TT 

6 A 2 

30 min. 

2 TT 

2 TT 

0. 8 TT 

0.8 tt 

V A1 







30 min. 

0. 1 T 

0.1 TT 

0.1 TT 

0.1 TT 


2 hr. 

12 TT 

0. 1 tt 

0.1 TT 

0 • 1 TT 

/ 

1 

34tt 

2 . Itt 

1. 3tt 

1. 2tt 

'J 6<j) 2 (T) 

/r 

sr 


/r 
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Figure 5. 
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